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SIKLSS INTENSITY FACTORS FOR 
BONDED ORTHOTROPIC STRIPS WITH CRACKS* 


by 

F. Dcldle and F. Erdoqan 
Lehigh University, Bethlehem, Pa. 

ABSTRACT 

In this paper the elastostatic problem for a nonhomogene- 
oiis plane which consists of two sets of periodically arranged 
dissimilar orthotropic strips is considered. It is assumed 
that the plane contains a series of collinear cracks perpendi- 
cular to Che interfaces and is loaded in tension away from and 
perpendicular to the cracks. First the problem of cracks fully 
imbedded into the homogeneous strips is considered. Then the 
singular behavior of Ihe stresses for two special crack geome- 
tries is studied in some detail. The first is the case of a 
broken laminate in which the crack tips touch the interfaces. 

The second is the case of cracks crossing the interfaces. An 
Interesting result found f » om the analysis of the latter which 
may have an important bearing on a possible delamination frac- 
ture initiation at stress-free boundaries in bonded orthotropic 
materials is that for certain orthotropic material combinations 
the stress state at the point of intersection of a crack and an 
interface may be bounded whereas in isotropic materials at this 
point stresses are always singular. A number of numerical 
examples are worked out in order to separate the primary mater- 
ial parameters influencing the stress intensity factors and the 
powers of stress singularity, and to determine the trends regard- 
ing the influence of the secondary parameters. Finally, some 
numerical results are given for the stress intensity factors in 
certain basic crack geometries and for typical material combina- 
tions. 

T 
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I. INTRODUCTION 

In considering the failure of a given structural component 
if the corresponding material is homogeneous and isotropic in 
its strength and thermomechanical properties, the related frac- 
ture process is relatively well-understood and the techniques 
dealing with such problems are sufficiently wel 1 -developed . 

This is particularly true in the absence of large scale plastic 
deformations around the dominant flaw from which the fracture 
failure would develop. On the other hand in composites, parti- 
cularly in fiber-reinforced laminates, the situation is much 
more complicated not only because of the nonhomogeneity and 
anisotropy of the material which make it very difficult to ana- 
1 yze the probl em, but also because of the highly nonhomogeneous 
and nonisotropic distribution of the strength parameter making 
the development and the application of a proper fracture cri- 
terion also very difficult. In such materials it is quite 
possible that the concept of the progressive growth of a domi- 
nant crack with a well-defined leading edge is not an appro- 
priate model for the cha rac ter i za t i on of gross fracture beha- 
vior. Very often the damage zone developing around the dominant 
flaw is somewhat irregular and diffused and the fracture process 
is generally governed by a principle of "weakest link", the 
local fracture propagation being progressive or in discrete 
steps. Nonetheless, whatever the gross mechanism governing the 
process of fracture failure in the structure, one may nearly 
always assume that locally fracture initiation and propagation 
will take place along the leading edges of the existing flaws 
where the conditions of the relevant fracture criterion are 
satisfied. Thus, in order to treat the local fracture pheno- 
menon in composite materials quantitatively, one may need the 
solution of the mechanics problem for flaws or cracks located 
at or near the phase boundaries or bimaterial interfaces. 

For composites which consist of bonded isotropic materials 
a wide variety of crack problems have been solved in which 
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t'ilht'f tht* .isvmptotic bohavior- of the stress slate arouml the 
points of neoiiietric sinqularity, or the results for a spec! tic 
I rack yeoiiietry tiave heen discussed (see, for examt)le, I and 
[2 \ tor nview and rt* f e rences ) . On the other hand, crack 
)rohlenis for homoneneous or nonhonio»ieneous anisotropic mat- 
erials remain relatively unexplored. Most ot the existing 
solutions refer to infinite planes 13-7). The crack problem 
for an orthotropic .trip is considered in [8] and that for an 
orthotropic sti'ip bonded to two orthotropic iialf plunes is dis 
cussed in 9'. The iletails of ttie problem for a finit(> crack 
located in the ne i qitborhood of, or i ntersect i no , a bimaterial 
interface in bonded anisotropic mater-ials do not seem to have 
been i nvi»s t i go ted . Cven though the problem is r'ather compli- 
cated mostly because of the large number of independent con- 
starrts entering the analysis, it may be made manageable under 
certain simplifying assumpt i ons . The mairr assumptions made in 
tins paper are (a) both materials are orthotropic, (b) the 
nonhomogeneous medium consists of two sets of periodicaMy 
arranged dissimilar strips having different thicknesses, and 
(c) the cracks in the strips are collinear, pe rpend i cu 1 a r to 
the interfaces, and also periodically arranged (Figure 1). 
Thus, one can take advantage of the symmetry of the medium 
and formulate the problem for two bonded strips only. The 

1. 0 r re spend i ng problem for isotropic layers or strips were con- 
sidered in *10] and (11]. In [12] the effect of Uio thickness 
and the elastic properties of the adhesive layer on the stress 
intensity factors i t> bonded dissimilar isotropic strips was 
cons i lie red . 

2. liFNERAl FORMULATION OF THE PROBLEM 

Consider the plane problem for an orthotropic medium. 
Referr'ing to, for example, [13] if u and v are the x and y com 
ponents of the d i s(> 1 a cement vector, the equations of equili- 
brium may be expressed as follows: 
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for plane strain. Here, E^j, , and G^j, (1,j) • (1,2,3), 
are the engineering elastic constants. Indices (1,2,3) refer 
to the (x,y,z) directions, and the matrix (b^j) given by 

(b,j) . B . A*' , A • (a,j) . (1 J) ■ (I .2.3) . 

'ii ■ ’'^i1 ' • <*> 

The s tress -d1 spl acement relations are 
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( 6 ) 


for generalized plane stress. 


Consider now the periodically arranged two sets of bonded 
orthotropic strips shown In Figure 1. In addition to the geo- 
metric symmetry Indicated in the figure, it will be assumed 
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th.it thr* mi'dium is suhjecti'ii to in-pl.in»» lo.ulimi which 
is syniini’tric with rospoct to th«» x axis ami is romoto 
from ami perpoiul i cu 1 a r to tho cracks. 

Thus, th»> solution of tho problem may be obtained by tne sian- 
dard superposition technique, from the viewpoint of fracture the 
important component beinq the perturbation solution in which the 
crack surface tractions are the only external loads. One may 
note that l)ecause of symmetry it is sufficient to consider the 
problem for one quarter of each strip only. Let (x^y) and 
(x 2 *y) be the local axes for tlie sets of strips 1 and ? as 
shown in Fiqure 1. Let the d i sp 1 acement s be expressed in terms 
of the following Fourier integrals: 


2 (" 2 r 

’ n] V" * iT] o-’j ^ 

Vj(Xj,y) » n!j(u,xp sin ya do ♦ -j nj(o,y) cos XjO do 


( 7 a , li ) 

where j » 1, and j » 2 refer to the strips 1 and 2, respectively. 
Substituting from (7) into (1) one obtains a system of ordinary 
differential equations for the unknown function, f^,..,nj, which 
are coupled in pairs. Solving these equations we find 


4 . . . ^ik ‘'’‘j 


f , ( o , X . ) • A . . ( o ) e 
J J k-1 






4 ^ 1 1 ''y / 'S K 

g.('^.y) * ^R.i(‘‘)o ^ , n,(o,y) 

J 1 J 

« ldj^R^^(o)e , (j*l,2) (8) 

In (8) Sji^ ,(j'l,2 , k*l,..,4) are the roots of the following 

c h a r a c t c r i s t i c equation: 
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4 2 2 

s * 


0 , s 


*J1 ’ "J4 


‘J2 • 


6j4 • • ^J5 “ ®J2^‘^J1 • ^ 

The functions Ajj^ and , (J-1,2 , k«l,..,4) are unknown 

and the constants Cjj^ and djj^ are given by 


‘^Jl " '*^j3 " • 


^i2 “ ‘^J4 


^’■^jl*j2^'®J3*j2 ’ 


‘^jl “ ■‘^J3 " ^*jr®Jl«j5^/^j3*Jl®j5 ’ 
^S2 “ “ ^*j2‘®J /j5^/^j3*J2®j5 • 


( 10 ) 


The unknown functions Ajj^ and which appear In (8) are 
determined from the boundary and the continuity conditions of 
the problem. In addition to the assumed nature of symmetry In 
loading and geometry, It should be emphasized that In the per- 
turbation problem under consideration the only external loads 
are the local self-equilibrating crack surface tractions. 
Consequently, both components of the displacement vector would 
vanish for y-*+«, and the x-component of the displacement, Uj 
(j*1.2) would be zero along the axis of symmetry Xj ■ 0 , 

(J*l,2). Thus, the sixteen conditions which have to be used 
to determine the unknown functions Ajj^ and Bjj^ (J"l,2; k*l,..,4) 
may be stated as follows. 

u^(hpy) = U2(-h2,y) , v^(h^,y) » V2(-h,,y) , 


0<y<« 


( 12 ) 
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a 2 xy("^’ 2 *y) * 0<y<» ,(13) 



Uj(0,y) • 0 , ° • 0<y<« , (j-1.2. , (14) 

Olyy(Xl .0) - p^ (x^ . |x^ |<a , 

v^(xpO) ■ 0 , a<|x^(<h^ , (16a, b) 

‘^2yy('‘2*0) " • 

V2(x2.0) ■ u , 0<|x2l<c , d<!x 2 l<h 2 (17a, b) 


In (9) It may arbitrarily be assumed that 
Re(Sj,)>0 , Re(s^2^'*0 • (J’1‘2) . (18) 

From (7), (8), (11) and (18) It therefore follows that 

Bji(a) - 0 , Bj2(a) • 0 , (J-1,2) . (19) 

Ten of the remaining twelve unknown functions may be eliminated 
by using the homogeneous conditions (12-15) In (8), (7) and (6). 
The last two unknown functions are then determined from the 
mixed boundary conditions (16) and (17). The problem may be 
reduced to a pair of Integral equations by defining 






0<|xj|<hj 


(J'1.2) 


( 20 ) 


and by replacing the conditions (16) and (17) by (20). Thus 
all the unknown functions Ajj^ and may easily be expressed 
in terms of the new unknown functions (^j and 4 » 2 * We now observe 
that part of the mixed conditions, namely (16b) and (17b) Is 
equivalent to 


(x^ ) * 0 


a<| x^ |<h^ 
0<|x2l<c 



( 21 ) 

0 .( 22 ) 


Substituting the results obtained from (6), (7), and (8) 
into the conditions (16a) and (17a) we obtain two integral 
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equations to determine and (^ 2 * 

Because of tne large number of elastic constants and un- 
known functions the process of deriving the Integral equations 
Is rather compi Icatod and 'engthy. However, the technique Is 
straightforward and Is quite similar to that followed In [10] 
and [11]. Therefore, the details of the derivations will not 
be given In this paper. As In [10], it can be shown that the 
Integral equations are singular and may be expressed as follows: 


+ t + x^^ (x^ ,-t)]^^ (t)dt 

1 


[ki2(X| ,t) 

L 


2 




- ki2^*l •-t)]<l>2(t)dt 

k 2 i ^^^2 ’"^^^*^1 ^ 

k22(^2*"^' * 


■ S7 


x,a, , 

‘•Tr't-x, t + x,^ 

(22 


x^cLo t (23a, b) 

C t. 

where and L 2 refer to the cracks on (y*0 , 0£x^<h^) and 

(y»0 , 0 <X 2 <h 2 ) In the strips 1 and 2, respectively, and 

“I ■ ^'^lyy''u/<'-“lxy''lyx> • “z ' ^^Zyy < ' ■'’Zxy''Zyx > ' 

^24) 

In deriving the Integral equations one needs to define In a 
systematic fashion a large number of elastic constants and 
Intermediate functions. Therefore, in order to conserve space 
the definitions leading to the expressions of the kernels 
k^ j , ( 1 ,j*l ,2) , and the constants and Yi^. and to the 
relationships between the functions Ajj^, and (}>^ will also 
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( » ) 

be omitted In this paper' These definitions and the details 

of certain derivations may be found In [14] for the group of 
orthotropic materials which would give a characteristic equa- 
tion having only real roots Sj|^,(J»l,2; k»l,..,4) (defined 
henceforth as the orthotropic materials of type I), and In 
[15] for materials which would give a characteristic equation 

with only complex conjugjte roots (defined as the orthotropic 

. I ** ) 

materl al s of type II)' ' . 


form : 


The kernels k^j which appear In (23) are of the following 


kij(xf.t) - k^j(x^.t.a) da , (l.j-1.2) . (25) 

Examining the behavior of K^j for a-*0 It can be shown that 

Kij • ^ ♦ 0(1) (26) 

where c^j are known constants. Even though this may Imply 
divergent kernels, by writing 

J '-J 

*[<,jdt [ -^-da . (l.j-1.2) (27) 

Jl j k 

and by using s 1 ngl e-va 1 ueaness conditions (see (21), (22)) 


f 

I 


U.(t)dt - 0 , (J-1,2) (28) 

‘■1 

It Is s'een that the singularity at a ■ 0 may easily be removed. 

Also, by examining the behavior of the Integrands K^j, (1,j»l,2) 

for a-»® It can be shown that they decay exponentially provided 

the series of colllnear cracks L. (1*1,2) are fully Imbedded In 
r*'] ' 

' 'Note that the constant Y14 is the same as the constants m,. 
and f,. defined in [8] (eqs.^16 and 19) and the constants 
defined In (24) correspond to 4u/(1+k) for the Isotropic materials. 
( * * 1 

' 'In practice, since p.. in the characteristic equation (9) 
appears to be always a negative quantity, the third type of mater- 
ial giving four pure Imaginary roots is not a realistic one. 
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the honioqeneous strips (i.e., they do not touch or intersect 
the bimaterial Interfaces). Thus. In solving the Integral 
equations (23). kfj may be treated as Fredholm kernels. In 
this problem since the kernels of the Integral equations have 
only a Cauchy type singularity, the functions would have j 
square root singularity at the end poIntsofL^and the equations 
may easily be solved by normalizing the Intervals and by using 
the technique described, for example. In [16]. After solving 
the integral equations, the stress Intensity factors may be 
obtained In terms of the functions For example, let Fig- 

ure 1 describe the crack geometry, I.e., let L]»(0,a), 
ic.d); then, the stress Intensity factors may be defined 
and obtained as follows [8]: 


k(a) • lim .TTt-TT o, (t.O) 
t*a 

k(c) ■ 11m /2('c“lT 
t^c 

k(d) ■ 11m /2 ( t- d ) (t.O) 

t-»d 

3. CRACK TOUCHING THE INTERFACE 


-lim Ml /T(a- tl 4»i(t) , 

t-a ' 

1 1m Mo/ 2( t-c) ^y(t) , 

t-c ^ ^ 

-11m t) <l> 2 (t) 

t-*d 

(29a-c) 


Two limiting cases of the problem discussed In the pre- 
vious section are physically Important and mathematically In- 
teresting. These are the cases of a broken laminate corre- 
sponding to a crack touching the Interface (e.g., a"h^, d<h 2 . 
Figure 1), and a crack Intersecting the Interface (e.g., a«h^ , 
d»h 2 , 0<c<h2» Figure 1). For example, referring to Figure 1, 
let a«h] and d<h2* In this case It may be shown that as 
and for -hi <( X] , t )<h] , c<( | X 2 | , 1 1 1 )<d the Integrands K 12 . K 21 , 
and <22 In (25) decay exponentially. Therefore, the kernels 
k^2* *^21* ^22 bounded In their respective closed do- 

mains. On the other hand for xi-^h], t-»h^ the exponential 
decay in Kii(xi,t,ci) disappears. Indicating that kii(xi,t) 
may contain terms which become singular as xi and t approach 
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the end point h^ s Imu I tiineous 1 y . These slnquler terms can be 
separated by studying the asymptotic behavior of the Integrals 
given by (25) (see [16] for the technique and [10], [14], [17] 
and [13] fo> the application). To give an Idea about the na- 
ture of these additional singular kernels let 

l^ll(^l«0 • * ^llf^*l*^^ * ^ ^ h^ , 

(30) 

where represents the singular terms and k^^ Is bounded 

in the related closed domain. Let the material be of type 1 
with the real roots (see equation 9) 


s ,, • .,>0 


0 


1 3 


-Ul 


1 


14 


Then the asymptotic analysis of (25) would give 

(h^-t)B^g/iD^'*’h^.i(^ 


k 1 ^ 5 ( X 1 , t ) • > 


85 
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+ \ 


88 
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15 /^l+hiu) 


“(oiXi )* 
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-t) 6 i g/uii 
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1 g/u)i •►hiu) 


- ( u )2 X 1 ) * 


(h, 

-t ) 6 i 5/<‘»2 

>hi 

'“I 

[(h,- 

•t)B 

1 g/ui 2 + hiu) 


“(u»iXi )* 


(hi 

-t)Bi 5 / 0)2 

*hi 

0*2 

[(hi- 

-t )6 

15/u.2*hiu) 

2 l‘ 

”(11)2X1 )^ 



0 < 

(x 

1 , t) < hi 


-Ui 


2 • 
(31) 


(32) 


where ^85’**'’^88 known constants and depend on the elas 
tic properties of the materials only [14]. 


i.Miether with I/(t-m), qlves a generalised Cauchy 

kernel. Substituting from (30) Into (23) and adopting the 
crack geometry shown In Figure 1 (with a«h]), the dominant 
part of (23) may be expressed as 


o 


-h, 

? ‘ ■ c<X2<d . 


-h^<x^<h^ , 


(33a. b) 


where the bounded functions and P 2 contain all the non> 
singular terms In (23). It Is clear that the solution of 
(33b) Is of the form 

■» 2 (t) - f 2 (t)[(t-c)(d-t)]'’^^ , c<t<d (34) 

giving the stress Intensity factors as defined In (29). The 
singular behavior of the solution of (33a) may be studied by 
1 p * < ng 


♦l(t) • F^(t)/(hJ-t*)^ . 0<Re(Y)<l . -h^<t<h^ (35) 


and by using the f unc t1 on - theore t1 c method described In, for 
example, [16]. Thus, If we define the following sectlonally 
holomorphic function 


1 fPl 'I'l ( t ) 

G(^) “ - J • (’‘l " <<e(2)) , (36) 

-h, 


by using (35) the asymptotic analysis of (36) gives 
F,(-h,)e^”Y , F,(h,) , 

G(i) . -i ! ! ! — !— ! 

(2hp^s1miY (z + h^^ (2hp^s1nnY (z-h^)^ 

1Gq(z)1< — ~ y- . Yo^Re(Y) . 


(37a, b) 
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where Cy and are real constants. Now, substituting from 
(37) and (32) Into (33a). we obtain the following characterls 
tic equation to determine the unknown constant y: 


-2cosny ♦ X 


I . 


- ZY 


li) 


85 ,-y 


♦ X 


1 


86 




(^>2 


* 87 - 


* • ( 58 ) 

-^15 

where, again the material type I Is assumed. It can be shown 
that for all material combinations (38) may have only one rooi" 

In the strip 0<Re(y)<l and this root Is always real. It can 
also be shown that as the orthotropic material constants tend 
to those for a pair of Isotropic materials, the root 7* obtained 
from (38) approaches the root of the corresponding Isotropic 
characteristic equation given, for example. In [10] or [18]. 

For this crack geometry the "stress Intensity factor" may 
be defined In terms of the cleavage stress 02yy neigh- 

boring material which, from the fracture viewpoint. Is the 
most Important stress component. To calculate this we note 
that (23b) gives the expression for O2yy(x2«0) for -b 2 iX 2 <.h 2 , 
that is in the uncracked as well as in the cracked portion 
of the '.trip. We also note that in the neighborhood of |x 2 |*h 2 
the singular behavior of 02yy governed by the density 

function and the singular part of the kernel k 2 ^ . As in ^ , 
it may be shown that for o'®, t-*h|, X 2 -*-h 2 . the exponential decay 
in K 2 ^ disappears, indicating that k 2 ^(x 2 ,t) miy contain terms 
which become singular as X 2 and t go to the end point simul- 
taneously. If we again let 

^21^^2’^^ * ^21s^*2’^^ * ^21f^^2*^^ * (89) 

the singular part of the kernel may be separated and may be 
expressed as 



. . 

■ ‘'0'[(h,.t)S,5/.,n.>,]>.(„,x,)> 

♦ X 

, 

[(hi”t)p^g/(»i2'*’rt2^*2^^”^^2^?^* 


where u^ Jnd ^re the positive roots dnd S 22 
characteristic equation (9) expressed for the strip 2 and the 
constants \ are defined in [14], Thus, for the purpose of ana 
lyziny the singularity Opyy "*‘*y expressed as 


'W«.' 

where P 2 Q contains all the nonsingular terms. Upon substitut- 
ing from (35) and (40) into (41), the asymptotic analysis gives 


k(h, ) 


(42) 


where 0 (x«) remains bounded as x«*-h- and the "stress intensity 
0 2 cl 


factor" k(h-|) is found to be 
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k(h,) ^ n*Hm (t ) 
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4. CHACk CROSS 1N(1 THl INTERIACl 


( 43 ) 


Conslil«*r MOW tho Ciisr of crack crossfno tho intorfaco. 

In this prohliMti tho intopral oi|uations (?3) arc still valiii 
with tho two iMul points of tho cuts and I , joiniiui at tho 
intortaco. lor oxainplo, j’oforrinii to riijiiro 1, lot a ■ hp 
d ■ h,,, and 0<c<h..ln this caso at tho ond point hj or 

x^ ^ hp all four kornols will havo singular torms. 

TIu* sinijular parts kp^ and k^p comintj from kp and k^^ aro 
soparat»*d and aro pi von hy (3?) and (40). Quito similar 
»*xprossions may oasily ho ohtainod for ^??s 

Tho dominant part ot tho systom of sinpular intopral opuations 
may thon ho oxprossod as 


C 

\ 


j" 








I I ■ ( - h I . h p , I 2 ■ ( c . h 2 ) • * i ' 7 ’ ^ ’ 1 . ? ) ( 4 4 ) 
whoro in tho analysis tho symmotry condition of 
-^'^("Xp) is usod. If wo now lot 
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^l(t) 


^2(0 


~ — r • 1 


Ov Re (^^ ,5 )< 1 , 


(45a. b) 


and ileMiit* the following sectionally holomorphic functions 


. G^(.) ■ . (4^a.h) 

asymptotic expressions for and may bo obtained as [16] 


1 Fl(-h,)e” 

G ( / ) . [ J 

(2h^ )^si niTp (z + h^ 


F,(h,) 
(2-h, ) 


* G^o(z) , 


G (2) . ] 1 

" (h2-c)®sinn(S {z-c)^ (h2*c)^sinn6 (z-h^)^ 


+ G 2 q(^) * (47a,b) 

where Gj^(J « 1,2) has a behavior similar to that of Gq(z) 
which is given by (37b). Noting that outside their respective 
cuts G^ and G 2 are holomorphic, substituting from (47) into 
(44), and following the procedure outlined, for example, in 
[IbJ (sec, also [14] for details) we obtain 

F2(c)cotn6“0 (48) 


- 0 . (i=l,2) 

1 ■ .i J J 


(49) 


where the coefficients In the functions f^j((0 depend on the 
elastic constants of the two strips only and are given in [14]. 


Since l' 2 (c) and f^j(^ij) 


(j»l,2) are nonzero constants, (48) 
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gives the known result 5 ■ 1/? ana from (49) we obtain the 
following c ha rac te r 1 s t 1 c equation to determine the power of 
singularity p: 

A(r) • " 0 • (1.J-1.2) , 0<Re(iO<l . (50) 

It Is also important to note that the end point values 
and r^(h 2 ) are not Independent and are related by 

where p Is the root of (50). An additional condition such as 
(51) is necessary to obtain a unique solution for the system 
of integral equations (23), since in this case there is only 
one s i ngl e- va 1 uedness condition which has to be satisfied by 
the displacement derivatives and i*> 2 « namely 

^ [\l^(t)dt + [ “ 0 • (52) 

^c -hi ^-h2 

A systematic study of (50) indicates that for all material 
combinations the characteristic equations may have either no 
root or only a single real root in the strip 0-Re(p)<l. Also, 

B » 0 is always a root and there are no other roots with Re(p) 

» 0. In the fcr*qoing analysis only the possib'lity of a power 

singularity is investigated. The results show that for certain 
material combinations (50) indeed has no root in 0<Re^6)<l, 
implying that for these materials at the intersection of the 
crack and the Interface the stress state would 'e bounded. 
However, this analysis does not prove that in such cases there 
may not be a weaker, namely a logarithmic singularity. To 

1 7 


investl'jdte this question In (45) we let ■ U and substitute 
the result Into (46). We would then obtain the followinq 
asymptotic relations: 

F,(h,) 

G^(z) - ^-log(z-hp - log(z + h^) + • 


G2(z) 


t2(c)c'"‘ , 


+ 


+ 


"Th2“c7 


1 og ( z -h 




(53) 


where and 62 ^ are bounded near and at the end points 

z “ i-lij and G^^ has a behavior similar to (37b) in the neigh- 
borhood of z » c. Substituting now from (53) into the integral 
equations (44) we obtain 

F2(c)cotit6»0, (54) 

log(h^ -X. )rg^ jFj(hj ) « R^(x^) , (i«l,2) (55) 

where and R^ are bounded functions and the constants g^j , 
(i.j*li2) depend on the elastic constants only. Equation (54) 
again gives the known result 6 * 1/2. For (55) to be valid 
at * (i* 1 . 2 ) the coefficient of singular terms must 

vanish, or we must have 

\ «i i “ 0 • (^"1 *2) . (56) 

Since Fj(hj) , (j»l,2) are nonzero constants, from (56) it 
follows that 
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I'jljl ■ “ 


(57) 


To show (57) analytically *;eenis to be Impossible. However, a 


systematic numerical analysis Indicates that for the material 
comb 1 nat ions having |& > 0 as the only acceptable power singu- 
larity (57) Is Indeed satisfied Identically. Furthermore, 
these studies also show that (56) always gives 


1 { h ] ) 1 ( ^ ) 

(hjT ■ Tfr‘^7 ■ •' 


The result expressed by (58) meaning that In the composite 
medium the derivative of the crack surface displacement is con- 


tinuous at the Interface Is, of course, the physically expected 


resul t . 


For the pair of materials In which (50) has a root In 
0 < 3 < 1 , at the point (y*0, x^»h^ or X 2 *-h 2 ) the stress state 
will be singular. At tills point, since the Important stress 
components are the normal and shear stresses on the interface, 
we may directly analyze the singular behavior of these stresses. 
To do this one has to go back to the original formulation of the 
problem and express these stresses In terms of the density 
functions and <P 2 . Thus, after somewhat lengthy but straight- 
forward analysis we find [14,15] 


1 




.y) 



(y .s)<|)j 


(s )ds 


(i*x,y) , 


= ( -h^ ,h^ ) , L 2 * (c ,h2) , (59) 
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studying the asymptotic behavior of the kernels h^^ It can be 
shown that as y*-0, s-* + h^ ln<-h^<s<h^, and s*h^ in c<s<hp simul- 
taneously the kernels become unbounded. By expressing 
h^j(y.s) » h^j^(y.s) ♦ » 

the singular parts h^j^ of these kernels can again be separated. 
For example, for obtain 

(h,+s)y,/2 (h,-s)Yi/2 

(h^+s)Y2'Y]]/2Y^2 (h^-s)Y2Y]^/2Y^2 

(h^+s)* + (u}2y/3^5)^ (h^ - s)^ + (u)2y/0i5)* 

+ ^3^81^’^4 ^82 (^ r^)Pis/‘^i 

^80 [(h^ - s)B^ 5/«*)1 ]*+y* 

+ ^3^83^^4^84 (^^rs)3is/u)g ^ 

^80 [(h^-s)6^g/a)2]*+y* 

where the definition of the material constants y and X as well 
as the expressions for the remaining functions h^^j* ^yis* 
hy 2 s "’^y found in [14] and [15]. 

If the materials are such that the stress state at 
(y*0, x^»hp is singular, i.e., 0<3<1 , then one can again de- 
fine and , (i»1,2) as in (45) and (46) and obtain (47). 

Now observing that outside the cuts and L^. specifically 
along the y-axis and G 2 are holomorphic, one can use (47) 
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to evaluate the singular terms In (59) (see, for example 
[16-18]). It can then be shown that 

o^xy^^py) • p * (61a, b) 

where the "stress Intensity factors" may be expressed In terms 
of the density functions as follows: 

''xx • “xx • 
t-*h^ 

•‘vw • 111" (hi-t)® 4 ),(t) . (62a, b) 

^ ^ t-*h^ ' 

The constants u^x i^xy known functions of the elastic 
constants and may be found In [14] and [15]. 

For the material combinations In which (j)^ and have 
no singularity at x^ » h^ , x^ • -h 2 , (i.e.. If B « 0 is the 
only acceptable root of (50)), since the kernels h^j have 
singular parts of the form (60), from (59) It Is not at all 
obvious that the stresses too would be bounded at the point 
(y»0, x^*h^. This question can be examined by substituting 
from (45) with B » 0, 6 ■ 1/2 Into (59) and by going through 
a routine asymptotic analysis, which yields 

"lxx<'’l-''> ■ y * C(y) 

aixy(hi .y) ' D(y) (63a, b) 
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where C(y) and D(y) are bounded functions. It turns out that 
In all material combinations for which p ■ 0, the constant 
®xx Identically zero; therefore, the stresses are bounded. 

Considering the fact that In Isotropic materials the 
stress state at the Intersection of an interface and a crack 
Is always singular (I.e., 6>0)» from the viewpoint of delamlna- 
tlon or debonding fracture the practical Importance of the pos- 
sibility of having bounded stresses at such locations In design- 
ing with certain orthotropic materials needs no elaboration. 

5. NUMERICAL SOLUTION 

In this paper the numerical results are obtained for 
several specific types of crack geometries. In the first 
group of solutions It Is assumed that the cracks are fully 
Imbedded In homogeneous strips and (see Figure 1) 

a<h^ , c “ 0 , d « ^*=^2 • 

The single crack, a*0, bf0ora/0b»0is considered 
as a special case. In this problem the integral equations (23) 
are solved by using the Gauss-Chebyshev Integration method [16] 
with L^ « (-a, a) , L 2 ■(-b,b) and under the si ngle-val uedness 
conditions (28). The stress Intensity factors are then obtained 
from (29a) and (29c) with d ■ b. 

In the second group of solutions It Is assumed that a » h^ 
and 0£b<h2. In this case the Gauss-Jacobi Integration method 
Is used to solve the Integral equations. The details of the 
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numerical method may be found in [16,|4or 18J. After obtain- 
ing the density functions the stress intensity factors are 
calculated from (29c) (with d*b) and (43). 

In the third group of solutions it is assumed that the 
crack crosses the interface, that is, a ■ hp d ■ h 2 , o<c<h 2 
(Figure 1). In this case for 0>O, the integral equations (23) 
are solved by substituting from (45) and by using the Gauss- 
Jacobi integration method. Here the additional conditions are 
(51) and (52). After obtaining and the stress intensity 
factors are determined from (29b) and (62) (see again [16], 
[14] or [18] for numerical procedure). 
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The elastic properties of the hidterials used in the nu> 
mcrical examples are shown in Table 1. Materials 3, 4, and 6 
are basically isotropic and the remaininq materials are ortho* 
Irupic. for the materials 1 through 8 the roots of the charac- 
teristic equation (9) are real, meaning that they are of type 
1. Materials 9 and 10 are of type II for which (9) has complex 
conjugate roots. The numerical results given in this paper are 
all for the case of plane stress. Table 2 shows the material 
combinations used in the numerical analysis. The table also 
shows the powers of singularity y and ti at the point of inter- 
section of the crack and the interface corresponding to a crack 
terminating at the interface (a*hp d<h 2 ), and that crossing the 
interface (a«hj, d^h^. oO), respectively (figure 1). Unlike 
the isotropic materials, the characteristic equations (38) and 
(bO) giving y and 8 in bonded orthotropic materials arc quite 
complicated. They contain six independent material parameters 
and hence do not lend themselves to a relatively simple syste- 
matic parametric study. However, once the material combination 
is specified / and p can be determined quite accurately. 


Even though it is very difficult to separate the material 
parameters which influence most of the values of p and y, and 
the stress intensity factors for the imbedded cracks, the cal- 
culations show that in this respect perhaps the most important 
single material parameter is the longitudinal stiffness ratio 
Eiyy/E 2 yy. In order to assess the effect of the remaininq ma- 
terial constants a rather large number of calculations were 
done by fixing E ] yy ^nd E 2 yy. by systematically varying one at 
a time the remaining six constants, and by calculating y, p, 
and the stress intensity factor k(a), the latter for Imbedded 
cracks in material 1 only. The general trend is as follows; 

As E^xx* ^Ixy* '*lxy medium 1 containing the crack) 

are increased, y and k(a) increase, and as E2xx’ ^2xy‘ '^2xy 
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Tabli* 1 Elastic 
used In 


No . 

E X x/ 
10'’N/m- 
( 10‘*ps1 ) 

Jyy/ 

(lO^psI) 

i 1(0) 

1 

55.16 

(8.0) 

1 70.65 
(24.75) 

2(0) 

134.45 

(19.5) 

31 .03 
(4.5) 

3(1) 

154.77 

(22.447) 

155.83 

(22.6) 

i 

4(1) 

167.55 

(24.3) 

1 70.55 
(24.75) 

' 5(0) 

. 

10.07 

(1.46) 

31 .03 
(4.5) 

: 6(1) 

• - J 

30. 34 
(4.4) 

31 .03 
(4.5) 

7(0) 

1 

44.82 

(6.5) 

1 55.14 
(22.5) 

8(0) 

34.48 

(5.0) 

6.895 

(1.0) 

9(0) 

21 .37 
(3.1) 

66.88 

(9.7) 

10(0) 

1 1 

17.24 

(2.5) 

17.24 

(2.5) 

— 


{*) The materials are boron-ep 
with various ply orientati 

25 


ts of the materials 
al calculations* 

















Table 2 The power of stress singularity y for a 
crack In medium 1 terminating at the 1n« 
terface and 0 for a crack crossing the 
Interface. The properties of materials 
used In various combinations are given 
In Table 1 


Comb . 

Material s 

Power of Sing. 

^lyy 

XV 

Med.l 

Med. 2 

y 

e 

Mm 

I 

1 

2 

0.55048 

0 

5.50 

0.20 

II 

3 

2 

0.65699 

0.04248 

5.02 

2.48 

III 

4 

2 

0.66549 

0.04887 

5.50 

2.58 

1 IV 

4 

6 

0.68914 

0.14547 

5.50 

5.75 

V 

4 

5 

0.80352 

0.05354 

5.50 

70.8 

VI 

7 

8 

0.74523 

0.05197 

22.5 

1 .40 

VII 

2 

1 

0.42258 

0 

0.182 

5.00 

VII I 

2 

3 

0.36911 

0.04248 

0.199 

0.403 

IX 

9 


0.61554 

0.08520 

3.88 

2.6 

j 

X 


9 

0.43410 

0.08520 

0.268 

0.384j 
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are Increased, y and k(a) decrease. Among these variables 
the most significant factor influencing y and k(a) appears 
to be the ratio of shear moduli ^ixy^^2xy' partly 

be observed also from Table 2 and Figure 2. The figure 
shows the stress intensity factor k(a) for Imbedded cracks 
in material 1 as a function of the width ratio h 2 /h^ for 
a fixed relative crack length a/h^-0.8 and for material 
combinations I, III, IV, and V given in Table 2. For these 


material pairs the stiffness ratio 
whereas Gixy/G2xy 2.58, 5.75 and 



70.8, 


constant 
respect! vely . 


It is seen that k(a) is consistently higher in material 


pairs having the greater ratio. Figure 2 also 

shows that for h 2 ** 0 , as expected, in all material combina- 
tions k(a) approaches the periodic collinear crack solution 


in an infinite plane which is the same for all homogeneous 
orthotropic as well as isotropic materials. 


A close examination of the results giving B. y. and k(a) 
indicates that generally one could accomplish a certain re- 
laxation in the stress singularity at the point of inter- 
section of a crack and an interface in composites by intro- 
ducing orthotropic materials. This may be seen, for example, 
by comparing the 6 values for various material combinations 
given in Table 2. In fact f(>: certain orthotropic material 
combinations it is even possible to have B*0, i.e., no 
singularity, whereas in isotropic materials 0 <b< 1 , i.e., the 
stress state is always singular. The value of B has, of course, 
an important bearing on the initiation of a possible delamination 
fracture from the stress-free boundaries in bonded materials. 

Even though the result regarding the possibility of 8*0 may 
appear to be somewhat paradoxial, considering the fact that 
in two isotropic wedges forming a half plane 8 is dependent 
on the wedge angles as well as the material constants and niay 
be zero for certain ranges of wedge angles, it should not be 
completel) unexpected. The possibility of reduction or com- 
plete elimination of singularity power 8 by varying the 
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secondary material constants seems to Introduce an added 
flexibility in designing against the edge delamination in 
bonded materials. 

In solving the integral equations it is assumed that the 
composite medium is under a state of generalized plane stress 
and is subjected to external loads away from and perpendicular 
to the cracks. Thus the crack surface tractions in the per- 
turbation problem considered in this paper are constant and are 
at the following ratio: 


PjtxT -l>2 


(65) 


The stress intensity factors obtained for the imbedded cracks 
located in the first or second set of strips are given in Fig- 
ures 3-7. Comparison of the results given in Figures 3 and 4 
shows that for the same longitudinal stiffness ratio Eiyy/E2yv 
and the same material 2, k(a) calculated for an isotropic me- 
dium 1 is consistently greater than that calculated for an ortho 
tropic material 1. This means that by introducing material 
orthotropy it is possible to obtain certain relaxation in the 
stress intensity factor. However, as seen from Figure 2, due 
to the effect of the secondary material parameters the opposite 
is also possible. In Figure 2 note that the combination IV 
refers to an isotropic-isotropic material pair - wher*. •. 3 III 
and V are i so trop i c -or tho tropi c pairs giving stress intensity 
factors which are respectively lower and higher than that of 
IV. Corresponding results for the stress intensity factor k(b) 
for cracks imbedded in the second medium are given in Figures 
6 and 7. Materials in Figures 5 and 7 are of the type II and 
those in Figures 2, 3, 4 and 6 are of type I. It siiould be ob- 
served that as the thickness of the uncracked strips go to zero, 
the stress Intensity factor in the cracked strips approach that 
of the periodic crack problem in the infinite homogeneous (iso- 
tropic or orthotropic) medium. 
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Piqui'C 8 shows a sample result for a composite medium in 
which both sets of strips contain cracks. Additional results 
for two as well as one set of cracks may be found in [14j and 

[15J. 

The stress intensity factors for the case of a broken 
laminate (i.e., for a*h^ and c"d, or a*0 and u»h 2 * c»0) are 
qiven in Fiqures 9-12. Figures 9 and 10 show the results for 
material combinations I and II where all materials are of type 
I and Figures 11 and 12 give an example for the material com- 
bination IX where both materials are of type II. The figures 
show that in all cases as the width of the uncracked strip 
(i.e., the net ligament between the cracks) goes to zero, as 
expected, the stress intensity factors become unbounded. In 
these problems the stress intensity factor is defined by (42) 
and is calculated from (43). 

The results for a crack crossing the interface are given 
in Figures 13-19. In these problems the stress intensity fac- 
tor at the crack tip k(c)»k^ is defined by and calculated from 
(29b). For those material combinations in which p>0 the stress 
intensity factors at he point of intersection of the crack and 
the interface kxx k^y are defined by (61) and are calculated 

from (62). For the material combinations II, IX and I used in 
these examples. Table 2 shows that power of stress singularity 
Y for a crack in mateiial 1 touching the interface is greater 
than 1/2. Therefore, as the crack length 2a, approaches 2h] or 
as c*h2, the stress intensity factor k^, at the crack tip calcu- 
lated on the basis of 1/2 power becomes unbounded. Also, as the 
length of the net ligament 2c goes to zero k^ again becomes un- 
bounded. These features of the solution may be observed from 
Figures 13, 16, and 19 giving the crack tip stress intensity 
factor as a function c/h 2 * Figures 13 and 16 show k^j for ma- 
terial combinations II and IX in which p>0. Figure 19 gives 
an example for the case in which p=0. It may be noted that 
qualitatively the results for the two cases are quite similar. 
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Iho stress intensity fiictors k^x *^xy material 
combinations II and IX are given in Figures 14, lb, 17, and 
18. Note that in the limiting case of c»h 2 . that is for the 
case of the crack touching the interface, the power of the 
stress singularity at the interface would be > which is al- 
ways greater than n. Therefore, as expected and as seen from 
the figures, for c »h 2 the stress intensity factors calculated 
on the basis of singularity power B become unbounded. In 
these problems for the type of loading under consideration 
the normal component of the stress intensity factor seems 
to be negative. Since there is no crack surface i n ter f erence , 
physically this means that normal stress along the interface 
near the crack surface is compressive, there is no inconsis- 
tency, and the singularity should be interpreted in the same 
way as in punch problems. 
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Figure 3 




ka/Pl 
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Figure 4 Same as Figure 3 for material combination II 




kb/P2'^ 



Figure 6 Stress Intensity fa 
in strip 2 (a»0, c« 
ation I 


J J 

h,/h2 2 


tor kb»k(b) for cracks imbedded 
, d“b<h2) for material combin- 
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Figure 9 The stress intensity factor k(hi)=k^ for broken 
laminates 1 in material combinations I and II 
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Figure 10 The stress Intensity factor k(h2)=kb for broken 
laminates 2 in material combinations I and II 




12 Same as Figure 10, material combination 






Figure 14 Normal component k^x of the stress intensity factor 
at the intersection of the cracks and the in- 
terfaces in material combination II 







Figure 16 Same as Figure 13 in material combination IX 




material combination IX 
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Figure 18 Same as Figure 15 in material combination I« 





O' 
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Figure 19 Same as Figure 13 in material combination 
where a=0 



